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Abstract 

The action of a nanoscopic spherically symmetric refractive index profile on a focused Gaussian 
beam may easily be envisaged as the action of a phase-modifying element, i.e. a lens: Rays travers- 
ing the inhomogeneous refractive index field n (r) collect an additional phase along their trajectory 
which advances or retards their phase with respect to the unperturbed ray. This lens-like action 
has long been understood as being the mechanism behind the signal of thin sample photothermal 
absorption measurements [Ullj, where a cylindrical symmetry and a different lengthscale is present. 
In photothermal single (nano-)particle microscopy, however, a complicated, though prediction-wise 
limited, electrodynamic (EM) scattering treatment was established [3j during the emergence of this 
new technique. Our recent study extended [4j this EM-approach into a full ab-initio model describ- 
ing the reality of the situation encountered and showed for the first time that the mechanism behind 
the signal, despite its nanoscopic origin, is also the lens-like action of the induced refractive index 
profile only hidden in the complicated guise of the theoretical Mie-like framework. The diffraction 
model proposed here yields succinct analytical expressions for the axial PT signal shape and mag- 
nitude and its angular distribution, all showing the clear lens-signature. It is further demonstrated, 
that the Gouy-phase of a Gaussian beam does not contribute to the relative photothermal signal 
in forward direction, a fact which is not easily evident from the more rigorous EM treatment. The 
model may thus be used to estimate the signal shape and magnitude in photothermal single particle 
microscopy. 
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I. INTRODUCTION 



Photothermal lens spectroscopy (PLS) has become a valuable tool in the study of solids 
and liquids [3H7]. Recent publications include the study of non-linear effects [8J and nanopar- 
ticles in solution [7]. Many authors have focused on the theoretical description of the thin- 
sample slab geometry which is often utilized in such macroscopic lensing experiments, often 
providing numerical equations which may be used to obtain absorption coefficients [1, 2j. In 
all these models the thermal lens induced originates from the absorbed power of a heating 
laser which constitutes a spatially extended cylindrically symmetric heat source in the heat 
equation. The solution obtained is then used to study the effect on the propagation of a 
probing laser beam. In thermal lens spectroscopy (TL) the probing beam is coaxial (possibly 
offset) with the heating beam [21 19] while in beam deflection spectroscopy [TO l ITT ] (BDS) the 
probing and the heating beams are aligned perpendicularly to each other. Both methods 
have in principle the same sensitivity jl2]. 

In contrast to these macroscopic techniques, a new microscopic approach has been devel- 
oped in the 1990's by Harada and Kitamori [T3| [Ti] termed photothermal lens microscopy. 
The developement of single particle photothermal microscopy [T514T7] followed and detects 
a very different kind of thermal lens with a mildly modified standard confocal fiuorescence 
microscope: Instead of an axially symmetric refractive index profile a spherically symmetric 
profile n (r), created by the point-like heat-source of the absorbing nano-particle, is probed. 
Also, instead of a profile that decays on the length-scale of the heating beam focus, in PT 
single particle microscopy a lens is probed which decays to half its value on the length-scale 
of the nanoparticle. On the other hand, the profile extends infinitely as 1/r and thus misses 
a characteristic length scale. This is the reason, why the description presented in this paper 
and within the recent EM-study [4j provide no evidence for a role-play of the Gouy-phase 
which is otherwise important for the probing of small scatterers as shown by Hwang and Mo- 
erner |T8j. While models for spherical absorbers have been put forward [T^ [20] . these were 
numerical in nature and targeted for large /xm-sized absorbers. The theoretical description 
of the nanoscopic photothermal lens has been first given by Berciaud et al. [3j in a scattering 
treatment relying on an extinction mechanism. Our recent ab-initio theoretical description 
of the electrodynamic problem has shown, however, that instead a simple lensing mechanism 
is responsible for the photothermal signal of nanoscopic absorbers showing a clear angular 
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thermal diffraction signature and a double-lobe lens signature in axial scans while disproving 
the assumption of interference-dominance in forward detection in the general situation. 

It is thus the aim of this paper to show that both worlds, the vast literature on macro- 
scopic thermal lens spectroscopy and the recent emerging tool of single particle thermal 
lens microscopy are very similar. The diffraction picture which has been a successful tool 
in the first domain will be shown to yield analytical and easily tractable expressions in 
the nano-scopic domain. They will allow for a quantitative assesment of absorption cross- 
sections of single nano-object based on standard photothermal measurements and will be 
able to explain the main phenomena of photothermal microscopy qualitatively as well as 
quantitatively, while providing an intuitive picture of the working mechanism. The quality 
of the simple model is checked against the more elaborate electromagnetic model within 
the extended scattering description. Axial scans and angular patterns of single heated nano 
particles will be described and compared. The angular diffraction pattern will be shown to 
explain the signal inversion observed for the first time for a single nano-particle upon the 
introduction of an inverse aperture in the detection path. 

II. THE LENS 

The lens to be considered in single particle generated nano-lens experiments such as pho- 
tothermal microscopy originates from the absorption of optical power provided by a focused 
laser beam by a single nano-particle. The absorbing particle can be treated as a point-like 
heat-source, yielding the steady-state temperature profile T(r) = To + AToR/r which decays 
with the inverse distance r from the particle of radius R. In the case of modulated heating, 
as utilized in the lock-in approach common to photothermal single particle microscopy, the 
assumption remains valid as long as the modulation frequencies used remain below 1 MHz, 
typically [31 H]. By the temperature dependence of the surrounding mediums' refractive 
index n, a corresponding refractive index profile n (r) is established: 



with An = ATq {dn/dT) being the heating induced refractive index contrast and no = 
n (To) the unperturbed refractive index. The amount of energy absorbed by the particle is 
determined by its absorption cross-section dabs and the intensity of the heating laser Ih (zp) 
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at the particle position. Together with the thermal conductivity k and the radius R of the 
particle this controls the induced temperature and thus the contrast An of the lens 

. (Tabs 4 / dn \ 
^"=W?(,dTj- 

The lens described by Eqn. [T] decays to half its maximum perturbation An at a distance 
of r = 2i?, making it a nanoscopic object. Nonetheless, it has infinite extent. It will only be 
limited by the probing beam which will for most cases be a confined and focused beam of 
diffraction limited extent. As An will quantify the photothermal single particle signal, the 
knowledge of the heating beam intensity and the thermal properties of the embedding bulk 
material will then allow the determination of the particles' absorption cross-section (Tabs- It 
is thus necessary to obtain expressions for the photothermal single particle signal from the 
induced refractive contrast An. This is the purpose of the next section. 



III. THE DIFFRACTION INTEGRAL 



The Frauenhofer diffraction integral for circlular apertures in the Fresnel-grade approxi- 
mation connects the complex field-amplitude of the probe beam Ua (p) in the aperture plane 
at z = to the complex field amplitude U (r, z) in the image plane at a distance z (see 
Figure |l| a)) [211 122] • A further fact or is included which represents the collected phase A$ 
for each component- wave as a result of the thermal lens (see Figure [l| b)) [23]. 



U (r, z) 



iz 



exp 



ikz 



Ua (p) exp 



R 



' 2z 



kpr 



exp (— iA$ (p)) pdp 



(3) 

To simplify matters, we will consider an axially symmetric scenario only, i.e. the refractive 
index profile will be symmetric with respect to the optical axis. This corresponds to the 
case of a heated nano-particle positioned along the optical axis (see Figure [l| a)). 

The phase advance A$ (p) for a single ray passing the lens immersed in a sample slab 
of thickness L at a distance p from the optical axis can be approximated via a straight ray 
calculation already utilized for coUimated beams in |23] : 
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Fig. 1. The scalar diffraction model, a) geometry for the diffraction integral (Eq. [3]). The shading 
in the image-plane corresponds to the difference in intensities between the diffraction of a cold and 
a hot nano-particle, Eq. jsj b) geometry for the computation of the phase advance (Eq. |4| 



where smh~^{z) = log (z + y/l + 2;^) and L ^ p were used after the integration. Although 
the integration in the distance p extends to infinity, the weighting by the Gaussian field 
amplitude Ua (p) will ensure the validity of the inequality. 

The field amplitude in the aperture plane will be taken to be the complex field of the 
probing Gaussian beam |24| having a focus displaced by Zp from the center of the heated 
particle (see Figure [T| a)): 
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where the beam- waist in the aperture plane is u (zp) = [l + Zp/z]^~^ , the curvature 



is given by Rc (zp) = Zp [l + z^/z^] and the Gouy-phase (g ' 



tan ^ {zp/ Zfj). The 



Rayleigh-range z^i is connected to the beam-waist Uq and the wave-number k = koUo by 
Zji = kuQ/2, where k^ = 27i/X represents the vacuum wave-vector for the wavelength A. 

To put photothermal single particle microscopy on a quantitative footing jl], we have 
introduced the relative PT signal S. It is the change in intensity / oc \Uf within the 
image-plane relative to the constant much larger background of the unperturbed field: 



S {r,z) = ^\U (r,z)|L(ATo(.,)) " 1^ (^'^)lL=o] / 1^^ ('^ = 0'^)lL=o (6) 
This relative signal will be independent of the probe power as long as no additional 



heating is induced, meaning that An = An (Ph) = const, with Pd- In case of a single laser 
being diffracted by its own induced thermal lens, the relative signal will be proportional to 
Pd- 

Through direct computation it can be shown that for small particles R < 100 nm the 
relative signal is not affected much by the geometrical diffraction of the inverse aperture 
disc of radius R. The diffraction-integral may be extended to i? = without an appreciable 
change in the result. We find the following analytic expression on the optical axis (r = 0), 
corresponding to the intensity- change detected in a closed aperture scenario: 

S (zp) = l^-^^A^'^or (1 + tRAnko)\^ - 1, (7) 

with the abbreviation ( {z,Zp) = (zp) + ik/ {2z) + ik/ {2Rc{zp))) and F (z) being 

the complex valued gamma function. The units of this expression seem to be odd at first, 
but a factor of the units meter [m] to the appropriate power has been omitted and ensures 
the unit less-ness of the expression. This expression being valid in the intermediate zone, the 
dependence on the image-plane distance z may be dropped in the far-field region by leaving 
out the second summand of (. The resulting dispersion-like signal shape (see Figure [2] a)) 
obtained through this expression may have been anticipated from the close similarity of the 
situation presented here and the thermal lens model used in thin sample slab geometries 
[H [21 [25]. In these cases, the parabolic refractive index profile approximation for example 
yields an optimal probing beam offset relative to the sample of one confocal distance, Zp = Zr. 
Here, the maximum signal is obtained when the heated particle is offset by about Zp = 0.7 Zr. 
For the off-axis image-plane signal, where the particle is still assumed to be positions on 
the optical axis, the radial coordinate r in the image-plane was transformed to an angular 
coordinate via tan (9) = r/z. The expression for the signal (Eqn. |6]) detected under an angle 
^ for r 7^ then reads: 
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iFi denotes the confluent hypergeometric function of the first kind. It may be related 
to the complex-order Laguerre polynomial via iFi (— az, 1, b) = Lia (6) for a G Reals. For 
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small detection angles 6' — )■ the expression reduces to the on-axis expression Eqn. [7] since 
tan (6) — 7- and iFi (ci, 1, 0) = 1 for any number Ci. For a probing beam positioned behind 
the lens {zp < 0), the angular pattern described by this equation shows a peak towards the 
center and an annular dip at larger angles (9c = 22°, see Figure |2] b)). For the case of a 
probing beam being positioned in front of the lens (zp > 0), the angular pattern changes its 
sign, relative to the previous scenario. Overall, the energy of the beam is only redistributed 
by the action of the lens, i.e. an integration of S [6, Zp) sin (6) from to 7r/2 gives zero. 
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Fig. 2. a) z-scan of the rel. photothermal signal for on-axis detection, i.e. = 0°, Eqn. [t] b) 
Angular pattern of the photothermal signal for Zp = —z^, i.e. Eqn.|8] The plot has been normalized 
to 1.0 on the optical axis. The parameters used in a) and b) the same as detailed in Figure |4j 

From this result we may, similar to the treatment given in |26], compute the total detected 
photothermal signal if a finite detection aperture is used. In this case it is not the modulation 
of the intensity detected on-axis or under a certain angle 6, but rather one needs to integrate 
the angular spectrum given in Eqn. M over the angular detection domain: 
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wherein the factor F compensates for taking the constant background \U (r = 0, -2)|^„=o 
in the integrand for the normalization in Eqn. |8] instead of the true Gaussian intensity. 
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A change of integration variables has been done from r to via tan (9) = r/z and the 
power Pr contained within a radius r at a distance z, which is taken large as compared to 



the Rayleigh-range zr, was used as given by Pr (z) = Pq[1 — exp (— 2r^/a;^ i^))]- Further, 
the expressions Pq = 7tIqUq/2 and I {z) = IqUq/u"^ (z) were used for the fraction of the 
power collected in the spherical cap area Agp^c as obtained by taking a constant on-axis 
intensity, Pd,io = ^spec-^(-2), and by using the correct varying intensity as determined by 
the Gaussian beam. Thus, if only the central bump of the diffraction pattern is collected, 
as it is usually the case when a dry objective collects the high-NA focused probe beam, a 
z-scan will show a clear change in sign. Indeed, a z-scan as described by Eqn. [7] and Eqn. 10 



shows a change in sign and a dispersion-like signal if An is changed according to the axial 
shifting of the absorbing particle with respect to the local heating laser. This is just what 
one would expect for the collected signal behind a (modifiable) diverging lens if probed by a 
focused beam and the lens were moved along the beam axis across the focus while its focal 
length were to change in such a way that it has its smallest |/| close to the probe beam focus 
and gets weaker, i.e. |/| larger, for large offsets. A further consequence of the given angular 
pattern of the photothermal signal is the possibility of an inversion of the photothermal 
signal. By using an annular aperture, only the annular region described by Eqn. [8] may be 
collected, and the total collected signal as described by Eqn. [TO] will be of opposite sign as 
compared to the use of a usual circular aperture (Objective NA). Indeed, the reversal of sign 
in the photothermal signal could both be observed experimentally for R = 60 nm AuNPs 



and confirmed in the calculations using the GLMT (see section IV). 

Equations ([T]) and (|8, 10) present the main results of this paper and give the background 
normalized intensity change S = AI/I detected either on the optical axis or detected under 
a finite angle 6 with respect to the optical axis upon the introduction of the lens n (r) (Eqn. 
[T|. Hereby, the particle/lens was assumed to be displaced by a distance Zp with respect 
to the probing beam- waist position (see Figure [T|. It seems worth mentioning that in this 
approach the Gouy-phase terms, while present in the fields in the aperture plane, cancel 
each other in the relative photothermal signal since they do not depend on the integration 
variable p and may thus be taken in front of the integrals. The found expressions were used 
to generate the plots in Figures |2] and |4] and allow the study of z-scans in photothermal 
microscopy setups similar to the thin sample slab studies in [HI [26H30] . To this end, the 
collection angle depends on the collecting microscope objectives' numerical aperture through 
6'jnax = arcsin (NAd/'n-m)? while the illuminating objective determines the beam waist(s). 
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IV. COMPARISON TO RIGOROUS VECTORIAL EM TREATMENT 



To compare the found results Eqn. ([7] and [8]) to a more rigorous solution of the problem, a 
full vectorial electromagnetic treatment will be used in the following. Therefore, the problem 
at hand may be expressed as the scattering of a shaped incident probe field interacting with a 
multilayered scatterer. The scatterer described by Eqn. ([T| may be viewed as an unbounded 
gradient refractive index lens (GRIN). The scattering of such an object has been studied in 
the literature and a publicly available C-code is attainable through reference |3lj providing 
the scattering coefficients of the GRIN a^v+i and 6iv+i when the scatterer discretized into 
concentric spherical shells (Figure [sj b)). The generalized Lorenz-Mie theory (GLMT) which 
is applicable to any spherically symmetric scatterer |32] positioned in an arbitrarily shaped 
beam [33H35| . We have therefore adopted and modified it to yield an analytical expression 
for the power Pd of the total electromagnetic field (E*, H*) detected in the far-field under 
a scattering-angle 6. This treatment has been verified by exhaustive comparison to single 
gold nanoparticle scattering and photothermal microscopy |4J. The quantity of interest is 
the projected Poynting-vector S* under a certain angle 6 with respect to the optical axis 
(Figure |3| a)): 

dPa {6, 0) = S* ■ dA = SidA. (11) 




Fig. 3. a) Schematic of the Poynting-vector integration within the GLMT framework, b) Dis- 
cretization of the refractive index profile n (r) into concentric spherical shells for the computation 

of the multilayer Mie scatter coefficients |31| . 
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Within the GLMT formahsm the total field is mathematically decomposed into se- 
ries expressions of the incidence field E* and H* with complex expansion coefficients g^, 
and an outgoing scattered field E** and H'^ described by complex scatter coefficients a„ 
and bn- Thereby, the total field consists of the incident field and the scattered field, 
V* = V* + V for V = E or H. The scattered far-field is commonly expressed via 
the scattering amplitudes through Eg = ^ exp (— ifcr) cos (0) 5*2 (^), Hg = —^E^, 
E'^= exp {-ikr) sin (0) Si {6), H'^ = §^EI. These amplitudes are: 



^ 2n + l 
2 ^-^ n (n + 1) 

n=\ ^ ' 



^ n„ (cos 61) _^ ^ (cos 61) 

" Tn (cos ff) " n„ (cos ff) 



(12) 



wherein the angular functions are 11™ (cos 0) = -P™ (cos 6') / sin 0, (cosO) = 
dP^ (cos 6) I dQ and are the associated Legendre polynomials. This artificially decom- 
poses the detected power into three parts which are not individually detectable in the scat- 
tering situation: dP^^ = dPinc + dPgca + dPoxt- Now, the time-averaged projected Poynting 
vector may be computed via (5*^) = ^Re{EgH^ — E^^Hg*), yielding three terms: 



2{Si) = Re {EIW; - ElHl*)+Re {E^Hl* - E;H^g*)+Re {EIH;* + E^^W; - E^H^g* - E^Hl*) . 

(13) 

The shape of the incidence field is determined by complex-valued expansion coefficients 
Qn, the so-called beam shape coefficients (BSCs). These coefficients reduce to (7„ = 1 for 
plane-wave illumination, i.e. regular Mie Theory. For a particle illuminated on-axis by 
a weakly focused Gaussian beam the modified local approximation (ML A, |36]) has been 
developed: 

9n (zp) = Q exp (-g {n - 1) (n + 2)) exp . (14) 

Herein, Q = {1 + is'-f)"^ with the beam-confinement factor s defined through s = 
ujq/{2zfi) and the defocussing parameter 7 = 2zp/ujQ describes the displacement of the 
particle relative to the beam- waist Wq, with Zp < Q being the situation where the focus 



is between the particle and the collecting objective. The result of Equation (11), with the 
integration over the azimuthal angle carried out, in the far-field can now be expressed as 
differential cross-sections in any forward/backward polar angle 9, dPcxt {d) = — dfToxt (^) -^0 
and dPsca (^) = ddsca (^) -^o, where the Gaussian beam focus intensity Jq = 2Pq/ttuJq was 
used: 
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da^c. (0)= ^{\Si (e) r + 1^2 (0) n , da,.t (0) = p [Re (M) Re (Su) + Im (M) Im (S,,)] . 

(15) 

The auxilliary functions 5*12 = 5*1 (0) + 5*2 (0) and M were introduced: 

oo Oil 

M{0) = y2 (cos e) + 7:„ (cos 0)] . (16) 

^ — ^ n{n + 1) 

n=l ^ ^ 

To obtain the total cross-sections, one needs to compute Usca = Jq dcTgca (0) sin (0) d^. 
The absorption may be computed from dabs = c^ext ~ CTsca with: 

C^sca = X] + 1^"'^ + ' = "p" 5Z + 1^"'^ + ■ 

n=l n=l 

(17) 

For a Gaussian beam (on-axis) one may calculate the integrated flux of the collected 
beam directly, i.e. / S* ■ dA with 25*^ = Re{ElH'^ — E'^Hq*)/2. The result may be written 
as a Cauchy-sum af^^ = ^ J2^=i '^tnc n with summands af^^ „ given by: 

^inc n = NmgmNn-m+ign-m+1 / [SmS„_m+l " (" 1)" AmA„_m+i] slu 6'd6', (18) 

m=l -^0 



wherein A„ = n„ — r„ and = n„ + r„. To ensure numerical stability for small angles, 
a direct recursive determination of A„ = n„ — r„ may be used [37J. The above equations 
may now be used to obtain the angular spectrum of the difference signal, i.e. the PT signal 
analogously defined to Eqn. |6l 



dph°t - dP^l'l [dc^sca {0) - dcTcxt (^)]"L+r^L+i " [ddsca {0) " ddext l-^JJa^.b^ 

= ^i^^;^^id = " Aa^ (9) 

mc inc V / 

To obtain the total relative PT signal collected with a finite detection aperture, the total 
cross-sections (Xsca.ext have to be used instead of the differential ones. This corresponds 



then to Eqn. 10 in the diffraction model. The predictions of both models are displayed in 
Figure g The best agreement is found in the low- focusing regime (wq ^ A). In the strong 
focusing case, the GLMT may be expected to deviate since the beam shape coefficients used 
here rely on the low focusing expansion of the coaxial field, while within the diffraction 
treatment the straight ray phase advance approximation (Eqn. |4]) becomes less applicable. 
Both frameworks agree for typical experimental parameters within a factor of the order of 
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unity. A direct scaling-free comparison of the on-axis diffraction formula Eqn. [7] and the 



Gaussian GLMT prediction with Eqn. 19 is shown in Figure |4]b), where the sohd thin black 
line corresponds to Eqn. [7] and the red dashed line corresponds to Eqn. [T9j The qualitative 
agreement is illustrated by further z-scan examples with finite numerical collection apertures 



see Figure H c), NA = 0.3 and NA = 0.75), i.e. Eqn. 10, scaled by 1.5 to match the 



GLMT predictions. The angular pattern in Figure |4| a) of the photothermal signal has two 
consequences: The relative change in intensity as compared to the unperturbed beam, i.e. 
the relative photothermal signal, is maximal, if the detection is on-axis or in a small angular 
detection domain around the forward direction (see Figure |4] b) and c) ). Further, both 
frameworks predict a signal inversion for the collection of an annular region of detection 
angles (Figure |4] d)), for example for [^^min, ^max] = [21°, 31°]. To test this prediction, a 
central beam-stop experiment has been carried out and is detailed in the following section. 
In summary, the comparison of the simple and intuitive diffraction-picture, as put forward in 
the preceding section, and the rigorous vectorial treatment of a Gaussian beam scattered by 
the GRIN lens shows a near-perfect qualitative agreement, while a quantitative agreement 
within a factor of order unity is found. 
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Fig. 4- Comparison of the diffraction (black) and Gaussian GLMT model (red). R = 10, ATq = 
100 K, no = 1.46, dn/dT = -3.6 x 10"^ A = 635 nm, ujq = 281 nm. A,, = 532 nm, wo,^ = 233 nm. 
The Diffraction model results have been scaled by a constant factor of 1.5 except for the d) where 
the factor is 0.8. a) PT signal angular distribution with positive (blue) and negative (red) signal, 
b) On-axis z-scan (NA = 0,9 = 0°) of PT signal. The black dotted curve is the unsealed prediction 
of the diffraction model, c) On-axis z-scan for a finite NA detection, d) On-axis z-scan with central 
beam stop (inverse aperture, see Figure |5]). The grey curve corresponds to no central beam-stop 

(NA = 0.75 from c) ). 
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V. SIGNAL INVERSION OF AN AXIAL SINGLE PARTICLE SCAN 



Having demonstrated the equivalence of the results obtainable within the generalized 
Lorenz-Mie framework and in the simple diffraction model, the following part of the paper 
will show a direct and obvious consequence of the above findings. Whereas the collection of 
the probe-beam in an angular domain around the forward-direction will show a dispersion- 
like lens-signature as displayed in Figure [2| the collection of an annular region can invert 
the detected signal. Indeed, the measurement of the photothermal signal of a single gold 
nanoparticle oi R = 30 nm radius shows this effect. To this end, a central beam-stop was 
introduced into the detection beam path. The collimated beam of diameter Do = 10 mm was 
thereby reduced to an annular ring of = 9 mm inner diameter corresponding to an angle 
of 6'min = arcsin (NA^Dj/ {Donm)) = 27° as given by the Abbe sine condition and the used 
numerical aperture of the detection objective. The illumination objective was a high NA oil 
immersion objective and the resulting probe and heating beam waists were UQ^d = 281 nm 
and Uo^h = 233 nm, respectively. The beams were offset in axial direction by Azj = 350 nm 
to achieve a symmetric signal configuration with the aberrated beams. The observed z- and 
xz-scans (Figure |6] a)) of the signal measured correspond very well to the calculated scans 
(Figure [6]b)) and clearly show the inversion of the signal. In addition to a pure change 
in sign of the signal, the zero-crossing shifts. Although these details are only reconcilable 
with the predictions of the exact GLMT description (see supplement of [4j), which includes 
the details of the aberrated beams, the basic principle and the prediction of the effect is 



understood based on the angular diffraction pattern described in section III 
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Fig. 5. Experimental configuration for tlie central beam-stop experiment used to measure the 
photothermal signal displayed in Figure |6j Apart from the central beam-stop, a detailed description 

of the experimental setup is given in reference [3]. 
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Fig. 6. Photothermal Single Particle Signal z-scan (left) and xz-scan (right) measured (a) and 
computed (b) with the exact beam shape coefficients in the GLMT framework for a full NA = 0.75 
detection (black curve on the left and upper pictures on the right) and with a inverse aperture up to 
an detection angle of 6^ (dashed blue on the left and lower pictures on the right). The parameters 
used in the calculation are the same as those in our reference |3]. 
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VI. CONCLUSION 



While a close relation of Gaussian beam diffraction and scattering has been shown by J. 
Lock et al. [3H] for the case of spherical dielectric particles, the more complicated scattering 
approach was thus far the only theoretical approach for the emerging technique of single 
particle photothermal microscopy. Although this ansatz may be used, the signature of a 
simple lensing mechanism (e.g. in z-scans) suggests a more intuitive model. Employing the 
scalar diffraction formalism common in thin sample slab absorption spectroscopy we have 
demonstrated that the signal obtained in photothermal single (nano-)particle microscopy 
can indeed be understood as the signal of a phase-modifying element, i.e. a lens, despite the 
microscopic origin of the latter. In contrast to the former, the spherical symmetry of the 
heated single-particle lens allows a direct analytical evaluation of the relative photothermal 
signal. The analytical model presented here could be verified in shape and absolute value 
when compared to experiments and a full electromagnetic vectorial treatment as published 
by the authors earlier. Further, it could be shown, that the Gouy-phase of a Gaussian 
detection beam does not contribute to the relative photothermal signal in forward direction. 
The understanding of the mechanism and angular distribution of the photothermal signal 
was then shown to explain the observed signal inversion upon the collection of an angular 
domain corresponding to the outer angles only. This introduces a simple and intuitive 
model for single particle absorption measurements and allows the quantitative assessment 
of nano-particle absorption cross sections. 
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